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Thi  ■ -■  ”•  


s report  examines  the  performance  of  the  QUICKE2  source  code.^thored  by 
T.A.  Dellin  and  C.J.  MacCallum  of  Sandla  LaboratorieSr^The  version  referred 
to  in  this  report  contains  an  IRT .modification  to  accommodate  photon  attenua- 
tion and  electron  emission  for  multiple  layers.  Numerous  comparisons  are  made 
of  QUICKE2  calculations  with  other  computer  codes  and  with  experiments  for  a 
variety  of  conditions.  Information  of  a practical  nature,  necessary  for 
modeling  problems  with  the  code,  is  provided.  Specific  input  requirements  are 
listed  and  outputs  are  interpreted  by  means  of  a sample  problem. 
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1.  INTRODUCTION 

QUICK1;2  (Ref  I)  is  an  analytical  computer  code  which  calculates  bulk 
and  vacuum  emission  currents  produced  when  materials  are  subjected  to  x-rays. 
This  code  was  developed  by  T.  A.  Dell  in  and  C.  J.  MacCallum  of  Sandia 
Laboratories  in  1974.  As  a refinement  of  earlier  QU1CKL2  versions,  it  has 
seen  extensive  use  in  the  radiation  effects  community  since  its  initial 
release.  This  report  is  intended  as  a descriptive  extension  to  Reference 
1.  It  reviews  in  some  detail  the  analytical  models  used  in  the  code,  input/ 
output  requirements,  and  compares  its  solutions  with  the  results  from  other 
codes  as  well  as  with  relevant  experimental  data.  The  version  examined  in 
this  document  contains  an  IRT  extension  to  the  original  QUICKE2  subroutines 
that  permits  the  application  of  QQrCKfi2  to  multiple-plate  geometries.  The 
calculations  are  carried  out  in  a single  computer  run  for  a given  spectrum 
and  plate  combination,  and  are  expedited  by  the  user-oriented  features  of 
simple  free-field  input  formats  and  versatile  output  plot  capabilities. 

These  and  other  code  features  are  described  in  detail  in  the  remainder  of 
this  report. 

Documentation  is  divided  into  sections  as  follows. 

Section  2 Description  of  the  Physics 

Comparison  of  QUICKE2  with  Other  Codes  and 
lixper  iment  s 

4 Description  of  the  Input 

5 Sample  Problem  and  Output  Description 

6 Computer  Requirements 

A summary  of  the  capabilities  and  limitations  of  QUICKH2  and  its 
multiple-plate  extension  is  shown. 

y % 

1.  T.  A.  Dellin  and  C.  .1.  MacCallum,  "QUICKh2:  A One-Dimensional 
Code  for  Calculating  Bulk  and  Vacuum  Emitted  Photo-Compton  Currents," 

Sandia  Laboratories,  SLL-74-0218,  April  1974. 
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Physical  Modeling; 

Geometry 

Materials  available 
Number  of  elements  per  plate 
Number  of  plates 
Photon  energy  range 
Photon  angle  of  incidence 

Typical  outputs 

Calculat ional  Considerations: 

Computer  time  requirements 

Number  of  energy  groups 
per  spectrum 


Infinite  plates 
25  elements,  1 ^ Z ^ 92 
<8 

No  limit 

1 keV  to  10  MeV 

0 to  180°  measured  from  sur- 
face normal  of  each  plate 

Photon  spectra,  electron  spec- 
tra, and  angular  distributions 
and  yields 

20  to  60  sec  7600  CPU  per  elec 
tron  emission  spectrum 

Photons:  50 
lilectrons:  34 
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2.  DESCRIPTION  or  Till:  PHYSICS 


A review  of  the  ]ihysics  employed  in  QUICKi;2  is  given  in  this  section. 
I.xtensive  descriptions  of  the  formulation  arc  given  by  Dell  in  and  MacCallum 
in  Reference  I. 

QUICKI.J  physics  can  be  divided  conveniently  into  the  categories  of 
photon  attenuation,  electron  creation,  and  electron  transport.  Highlights 
of  each  segment  are  listed  below. 

2.1  PHOTON  ATTENUATION 

QU1CKE2  calculates  the  attenuation  of  photon  beams  passing  through 
materials,  assuming  that  the  flux  decreases  exponentially  at  the  distance 
penetrated  for  each  energy  bin.  Compton-scattered  photons  and  other  secon- 
dary photons  are  neglected.  The  energy-dependent  attenuation  lengths  are 
taken  from  the  tables  of  Biggs  and  l.ighthill  (Ref  2). 

The  calculation  proceeds  as  follows,  A photon  spectrum  is  input  into 
QlJICki;2  as  a set  of  coordinates  (JN/dE  or  dli/dE,  EJ . In  addition,  the 
composition  and  thickness  of  one  or  more  attenuating  layers  is  specified. 
For  materials  of  one  element,  an  attenuation  length  is  interpolated 
from  the  tables  for  each  idioton  energy  E.  The  number  of  photons  of  energy 
E that  leave  the  material  is  then  calculated  as  the  incident  number  times 
the  attenuation  factor  |exp(-T/Aj,) ] , where  t is  the  thickness  of  the  mate- 
rial. For  materials  of  two  or  more  elements,  the  i)hotons  are  assumed  to 
pass  through  a scries  of  monoelemental  layers  of  thickness  where  f. 

is  the  fraction  by  weight  of  the  i—  element  of  the  compound  and  t is  the 

2 

thickness  of  the  compound  in  g/cm  . The  number  of  photons  of  energy  E 
leaving  the  material  is  then  taken  as  the  incident  number  times  the  atten- 
uation factor  |exp(-Zjf j/Xj  . )t  1 . 


2.  (I.  E.  Simmons  aiid  J.  11.  Ilubbell,  "Comiiarison  of  Photo- Interaction 
Data  Set  VI 1.  Biggs-Eighthi 1 1 (Rev.)  and  ENDF/B,"  NBSIS  73-241,  duly  1973. 


Except  near  the  photoelectric  edges,  log-log  interpolation  of  the 
cross  sections  taken  from  the  tables  is  accurate  since  the  intervals 
between  photoelectric  edge  energies  are  closely  approximated  by  line  seg- 
ments for  a log-log  plot  of  cross-section  versus  photon  energy.  If  an 
interval  between  two  tabulated  points  contains  both  a spectral  energy  E and 
a photoelectric  edge,  the  relationship  of  E to  the  edge  energy  is  determined 
and  the  cross-section  is  interpolated  accordingly. 

The  tabulated  cross-sections  include  photoelectric  absorption, 
coherent  and  incoherent  scattering,  and  pair  production.  All  of  these 
processes  are  considered  to  remove  photons  from  the  incident  beam.  Eor 
photon  energies  relevant  to  QUICKE2  (1  keV  to  10  MeV),  the  error  intro- 
duced by  including  coherent  scattering  in  the  absorption  cross  section 
is  small  since  this  process  is  dominated  by  the  photoelectric  effect  at 
low  energies  and  by  incoherent  ( *■  Compton)  scattering  at  high  energies. 

2.2  ELECTRON  CREATION  AND  TRANSPORT 

Electron  creation  in  QUICKE2  is  due  to  Compton,  photoelectric,  and 
Auger  processes  occurring  more  than  an  extrapolated  electron  range  away 
from  any  material  interface.  Photon  energies  from  I keV  to  10  MeV  can  be 
specified.  Electron  emission  due  to  a spectrum  of  photon  energies  is 
obtained  by  superposition  of  emissions  from  monoenergetic  photons. 

Electron  angle  and  energy  distrubtion  functions  for  the  bulk  material 
are  obtained  from  exact  solutions  to  the  transport  equation.  They  are 
exact  only  for  the  regions  of  the  target  that  are  further  than  an  extra- 
polated electron  range  from  the  material  interfaces.  Distribution  func- 
tions for  electrons  emerging  from  infinite  slabs  of  materials  into  vacuum 
are  obtained  from  these  exact  bulk  solutions  using  interpolative  approxi- 
mations. Complete  emission  electron  energy  and  angular  distributions  can 
be  obtained  for  normal  photon  incidence,  and  for  non-normal  photon  sources 
with  azimuthal  symmetry.  Angle-integrated  spectra  are  available  for  other 
non-normal  photon  incidences. 
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COMI’ARISON  01-  QUlCKi;2  WITH  OTIil.R  COOLS  AM)  L.XI’LRIMLNT 


QIIICKI,2  calculations  of  electron  emission  are  compared  with  those  of 
other  codes  and  with  experiment  in  this  section.  Specific  comparisons  are 
made  to  the  SANOYL,  I’OLM,  (;i0\l'2,  and  LI  lUN  codes  (Refs  3-6)  and  to  the 
experiments  of  Barlett  and  Weingart  (Ref  7),  Bradford  (Ref  8),  Denison  et 

Ref  4),  and  Lbert  and  Lauzon  (Ref  9).  I’lioton  attenuation  results 

ompa  red  to  those  of  the  l-TIULL  code  (Ref  10). 
ihle  3-1  contains  total  reverse  yields  from  AH,  Cu,  and  Ta  for  a 
'-keV  brensstrahlung  spectrum  as  measured  by  Bradford  and  as  calculated 
by  the  POIiM,  S;\NI)YL,  and  QUICKL2  codes.  This  table  is  adapted  from  Ref- 
erence 7. 

Table  3-2  contains  forward  and  reverse  yields  for  5-,  10- , and  15- 
keV  blackbody  spectra  incident  on  40  mils  of  AH  as  calculated  by  S/\NDYL 
and  QUICKL2. 

Figures  3-1  through  3-3  show  energy  distributions  of  the  back- 
emitted  electrons  for  the  Bradford  experiments  and  for  QUICKL2. 

3.  H.  M.  Colbert,  ".SANOYL:  A Computer  Program  for  Calculating  Com- 
bined Photon-Llectron  Transport  in  Complex  Systems,"  Sandia  Laboratories, 
SLL-74-0012  (1974). 

4.  W.  L.  Chadsey,  "A  Monte  Carlo  Photo-Current/Photo-Lmission  Com- 
puter Program,"  National  Symposium  on  Natural  and  Man-Made  Radiation  in 
Space,  Las  Vegas,  March  1971. 

5.  li.  P.  Wenaas  and  M.  Williams,  "GRAPUP:  A Computer  Code  for  Sec- 
ondary Llectron  Lmission  from  Plates,"  Gulf  General  Atomic,  GAMO-10155, 
.June  1970. 

6.  "LTRAN:  Monte  Carlo  Code  System  for  Llectron  and  Photon  Transport 
Through  Extended  Media,"  NBS  reports  9836  and  9837. 

7.  C.  .J.  MacCallum  and  T.  A.  Oellin,  .1.  Appl.  Phys . 44,  1878  (1973). 

8.  .J.  N.  Bradford,  "Radiat ion- Induced  Llectron  Emission,"  AI'CRL-TR- 
74-0583,  November  1974. 

9.  P.  Ebert  and  A.  Lauzon,  "Measurement  of  Gamma-Ray- Induced  Sec- 
ondary Electron  Current  from  Various  Elements,"  ILEL  Trans.  Nucl.  Sci. 

NS- 1 3 , p.  735  (1966). 

10.  M,  .1.  Nowak,  M,  R.  White,  and  .1.  P.  Wondra,  "The  LIllIlF  Code 
for  X-Ray  LYiergy  Deposition  with  fluorescent  Heating  and  fluorescent 
fluxes,  Gulf  General  Atomic,  GAMO-9463,  .July  1969. 
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Figures  .'i-4  through  3-6  give  energy  distributions  for  reverse  emis- 
sion from  A?,,  Ti  , and  I'e  as  measured  by  Penison  et  al  . and  as  calculated 
by  POIiM  and  QU1(;KF2.  Because  the  absolute  yields  for  the  exiieriment  and 
I’OFM  were  not  available,  the  curves  were  normalized  to  the  same  electron 
number  densities  as  the  QUICKi;2  results  at  the  peaks  of  either  the  K-photo 
edge  or  L-Auger  contribution,  whichever  was  greater.  QU1CKF,2  yields  are 
given  on  the  figures,  and  various  contributions  to  the  electron  spectra  are 
labeled.  These  figures  are  adai)ted  from  Reference  4. 

Figures  .3-7  through  .3-10  show  the  reverse  and  forward  cumulative  emis- 
sion s|ioctra  from  typical  satellite  materials  for  2-,  5-,  10-,  and  15-keV 
blackbody  photon  spectra  as  calculated  by  t'iRAP2  and  QUICK1:2.  The  reverse 
emission  is  from  a quartz  solar  cell  cover,  while  the  forward  emission  is 
into  an  aluminum  equipment  box.  The  incident  spectra  used  in  the  latter 
calculations  arc  the  blackbody  spectra  filtered  through  a typical  solar 
cell  and  through  the  wall  of  the  box,  as  shown  in  Figure  .3-11.  In  these 
instances,  the  ordinate  label  "per  calorie"  always  means  "per  calorie 
incident  on  the  Front  face  of  the  solar  cell." 

In  Figures  5-12  through  3-15,  [ihoton  attenuation  predicted  by  the 
subroutine  ATTFN  is  compared  to  ITIUFF  calculations  for  the  cases  of  Fig- 
ures 3-7  through  3-10.  Here,  the  photon  spectra  are  also  comiiarcd  as  they 
leave  the  solar  cell  before  entering  the  aluminum  box.  The  quantity  jUot- 
ted  is  the  cumulative  fluence  of  photons  with  energy  less  than  the  energy 
indicated  on  the  abscissa.  Very  minor  differences  arc  seen  in  the  results. 

For  a hard-spectrum  comparison,  the  forward  and  reverse  yields  and 
the  forward  angular  distributions  calculated  liy  QU1CKF.2  arc  compared  to 
the  experimental  results  of  F.bert  and  Lauzon  for  1.25-MeV  photons  inci- 
dent on  thick  targets  of  C,  AH,  Cu,  and  Pb . The  results  are  summarized  in 
Table  3-3  and  in  Figures  3-lb  through  3-19.  Note  that,  while  the  forward 
yields  obtained  by  the  two  methods  agree  well,  the  QUICKF2  angular  distri- 
butions are  not  as  sharply  peaked  at  0 = 0 as  the  experimental  ones.  Fig- 
ure 3-20,  taken  from  Reference  7,  shows  Dell  in  and  MacCallum's  comparison 
of  QUICKr;2  with  the  F.bert  - Lauzon  experiment  and  with  the  Monte  Carlo  code 
SANDYL.  Their  results  are  in  units  of  electrons  per  unit  angle  rather 
than  electrons  per  unit  solid  angles,  emphasizing  the  fact  that  the  solid 
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;inj;le  approaches  zero  as  0 approaches  zero.  Since  this  form  is  indicative 
of  real  spatial  emission,  the  discrepancies  in  the  results  at  small  angles 
shown  in  figures  3-l(i  through  .3-1*.)  probably  have  minor  impact  on  calcula- 
tions employing  QUlCKfiZ  electron  emission  data.  The  results  in  Table  3-1 
through  3-3  and  figures  3-1  through  3-b  indicate  that  the  QlJlCKhJ  calcula- 
tions agree  quite  well  with  SANDVl.  and  I’OliM  Monte  Carlo  calculations  and 
with  the  experimental  data  of  Bradford.  The  largest  differences  were  for 
reverse  emission  due  to  low  energy  photons  ('Table  3-2 J but  the  discrepancies 
were  less  than  a factor  of  two. 
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Table  3-1 


QlJICKi;2  RT.SIILTS  COMPARIT)  WITH  liXl’l-.R IMHNT  AM)  OTIlliR 
CALCULATIONS  TOR  RKVT.RSi;  YIT.LDS  IROM  1)  I TFliRTNT  MATTRIAI.S^ 
(Photons  from  50-kcV  brcmsst rahlnng) 


M 

Cij 

la 

Bradford  [experiment 

1.9 

X 

1(  ■ ' 

4.8 

X 

10-' 

1 .2 

X 

H)-"^ 

(ATCRL) 

QUICKT2 

1 ..S 

X 

10’' 

3.0 

X 

10-' 

1.0 

X 

10-^’ 

-7 

-7 

•7 

SANDYL 

1.5 

X 

10 

3.2 

X 

10 

8.7 

X 

10  ' 

POT.M 

1.5 

X 

10'' 

3.3 

X 

10-  ' 

8.8 

X 

H)-' 

'*Trom  Ref.  7. 


Table  3-2 


QUlCKt:2  RTSULTS  TOR  TORWARI)  AND  RTVTRSi:  YIT.LDS 
FROM  M COMPARTD  WITH  SANDYL  TOR  DITTTRKNT  SPT.CTRA 
(Plate  thickness  40  mils) 


5-keV  BB 

10-keV  BB 

15-kcV  BB 

Reverse 

-7 

-8 

QUICKT2 

3.4  X 10 

9.2  X 10 

3.9  X 10 

-7 

-8 

SANDYL 

2.3  X 10 

6.6  X 10 

2.8  X 10 

To  rwa  rd 

-8 

-8 

, -8 

QU1CKT.2 

2.9  X 10 

2.9  X 10 

2.0  X 10 

-8 

-8 

. ^ .,,-8 

SANDYL 

3.4  X 10 

3.1  X 10 

1 . / X 10 

5-,  10-,  and  15-keV  blackbody  spectra  are  incident. 
All  yields  are  in  C/cal  incident  on  front  face  of  AX. 


Incident 

Front 
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Tabic  3-3 

QlIICKi;:  ,V\I)  liXPliRIMIA'TAI.  RbSULTS  COMPARi;!)  FOR  lllCIl-IiNFRGY  PII010NS 


Thickness 

F.xperimcnt  of 

Materia  1 

(K/cm-  ) 

F.bcrt 

1. 

Lauzon 

qi)ickf;2 

Forward 

Yields^ 

Carbon 

0.157 

3.0 

X 

10-« 

2.8  X 

io-« 

{P.  1 

X 

10‘'^1 

(8.4  X 

lo'-^ 

Alum i num 

().bl4 

0 *7 

X 

10-« 

2.3  X 

!()-« 

(8.0 

X 

10-^ 

(7.0  X 

Copper 

0.4b7 

1.8 

X 

io-« 

1 .8  X 

io-« 

(5.3 

X 

lo'-^ 

(5.5  X 

10-^ 

Lead 

0.591 

2.0 

X 

10"^ 

1.7  X 

10'^ 

(b.l 

X 

lO*'') 

(5.3  X 

lo-'b 

Reverse 

Yie Ids^ 

Carbon 

1 . 1 

X 

10-^-^ 

1 .2  X 

10-^ 

(3.4 

X 

10“^  1 

(3,b  X 

“4 

10  ) 

Aluminum 

1 . 1 

X 

10--’ 

2.2  X 

10'‘-^ 

(3.3 

X 

10'^ 

( fa  . ()  X 

10’^) 

Copper 

4.3 

X 

10-^ 

3.7  X 

lo'^ 

(1.3 

X 

lo'-'*') 

(1.1  X 

10'^) 

Lead 

7.3 

X 

-9 

10 

7,4  X 

10-^ 

(2.2 

X 

lO"^) 

(2.2  X 

10"^) 

‘S'iclds  without  parentheses  are  in  C/cal ; yields  in 
parentheses  are  in  elcctrons/photon . 
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(electrons/cdl/MeV) 


ELECTRON  ENERGY  (keV) 


igure  3-1.  QIJICKF.2  results  compared  with  experimental  results 
of  Bradford  for  reverse-emitted  electron  energy 
spectra  obtained  from  5()-kcV  bremsst  rah  lung  on  A*. 
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(electrop.s/cal/MeV) 


(electrons/cal/MeV) 


RT 


Figure 


ELECTRON  ENERGY  (keV) 


3-4.  QUICKFZ  reverse-emitted  electron  energy  distribution 
conpared  with  POEM  code  and  experimental  results  of 
Denison  et  al.  for  8-keV  monoenergetic  photons  inci- 
dent on  AH;  curves  normalized  to  same  heights  at 
K-photo  edge  peaks 
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Figure  3-5.  QUICKE2  reverse-emitted  electron  energy  distribution 
compared  with  POEM  code  and  experimental  results  of 
Denison  et  al.  for  8-keV  monoenergetic  photons  inci- 
dent on  Ti;  curves  normalized  to  same  heights  at 
K-photo  edge  peaks 
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(electrons/cal/MeV) 


Figure 


ELECTRON  ENERGY  (keV) 


3-6.  QUICKE2  reverse-emitted  electron  energy  distribution 
compared  with  POIiM  code  and  experimental  results  of 
Oenison  et  al . for  8-keV  monoenerget ic  jihotons  inci- 
dent on  Fe;  curves  normalized  to  same  heights  at 
L-Auger  peaks 


number  of  electrons/cal  with  energy 


Figure 


RT-13690 
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ELECTRON  ENERGY  (keV) 


100 


3-7.  Comparison  of  QUICKli2  and  GRAP2  electron  spectra 
for  reverse  emission  frcm  a solar  cell  cover  and 
forward  emission  from  AS.  after  attenuation  by  the 
solar  cell;  incident  photon  spectrum  due  to  black- 
body  radiator  with  temperature  2 keV 
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^^QUICKEZ 


Figure 


3-9.  Comparison  of  QUICKF2  and  GRAP2  electron  spectra 
for  reverse  emission  from  a solar  cell  cover  and 
forward  emission  from  M after  attenuation  by  the 
solar  cell;  incident  photon  spectrum  due  to  black 
body  radiator  with  temperature  10  keV 


number  of  el ectrons/cal  emitted  with  energy 


ELECTRON  ENERGY  (keV) 


Figure  3-10.  Comparison  of  QUICKF;2  and  CIKAP2  electron  spectra 
for  reverse  emission  from  a solar  cell  cover  and 
forward  emission  from  M after  attenuation  by  the 
solar  cell;  incident  photon  spectrum  due  to  black- 
body  radiator  with  temperature  15  keV 
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RLVTRSE  EMlSSinu 


4.0  X 10'^  g/cm^ 

6.6  X 10'^  g/cm^ 
2.4  X 10'^  g/cm^ 
4.9  X 10  ^ g/cm^ 

2.6  X 10'^  g/cm^ 


EQUIPMENT 

BOX 

CAVITY 


RT-06431 


Figure  3-11.  Materials  used  in  emission  and  attenuation  cal- 
culations to  represent  typical  satellite  solar 
cell  and  equipment  box  cover.  The  QU1CKE2  ver- 
sion described  here  performs  the  calculation 
with  a single  computer  run. 
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Figure  3-15.  Comparison  of  QUlCKIi2  and  FHIH'I-  15-keV  blackbody 
photon  spectra  attenuated  by  solar  cell  and  30 
mils  of  M 
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Figure  3-16.  QUICKF2  and  experimental  results  of  l-bert  and 
I.auzon  compared  for  forward-emi ssion  electron 
angular  distributions  due  to  1.25-MeV  photons 
incident  on  C 
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I ' 
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RT-13700 


20  40  61 

EMISSION 


Hisure  3-18.  QIJICKF.Z  and  ex])eriment; 

l.auzon  compared  for  foi 
angular  distributions  i 
incident  on  Cu 
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Figure  3-19.  QUICKF2  and  experimental  results  of  Fbert  and 
Lauzon  compared  for  forward-emission  electron 
angular  distributions  due  to  1.25-MeV  photons 
incident  on  Pb 
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4.  DHSCkIPTlON  01-  nil;  INPUTS 

Uctailcd  descr  ijit  ions  of  the  input  cards  rccjuired  h>'  QIJICKi;2  arc 
given  in  this  section.  \’ariablc  names  and  physical  or  calcu  lat  ional  sig- 
nificance are  given,  as  well  as  the  fomiat  for  reading  them  into  the 
code.  Most  of  the  variables  are  in  NAMHLIST  format.  These  are  specified 
in  free  form,  separated  by  commas.  The  only  restrictions  are  to  begin 
the  list  with  bSN/\MI.,  leave  column  1 blank  on  all  cards,  and  end  the  list 
with  bSbP.NU.  "b"  indicates  a blank  space  (necessary  only  in  column  1) 
and  "NAMP."  is  the  name  of  the  NAMP.LIST  given  in  the  variable  descriptions 
(sec  Figure  5-1  for  a sample  of  the  input  deck). 

In  the  input  descriptions,  arrays  are  indicated  by  an  index  following 
the  variable  name.  If  no  index  appears,  only  a single  value  is  read  in 
for  the  variable.  The  code  counts  the  numbers  of  variables  in  cases  where 
arrays  are  involved,  so  the  number  of  values  being  read  in  need  not  be 
specified.  Default  values  and  maximum  or  minimum  numbers  of  values  are 
given  where  appropriate.  blements  available  in  ()UlCKh2  are  listed  in 
Table  4-1.  Elements  not  listed  can  be  referenced. 

The  10PT(27)  option  described  in  the  table  causes  a file  of  data  to 
be  written  which  can  be  treated  directly  by  the  PLOTAl.l.  code  (Kef.  11). 
This  permits  the  overlay  of  photon  and  electron  spectra  from  different 
plates  or  from  different  QU1CKF.2  calculations. 


11.  A.  .1.  Woods,  T.  N.  Delmer,  and  M.  A.  Chipman,  "The  Arbitrary 
Body  of  Revolution  (ABORC)  Code  for  .SCCMP/ 1 IMP , " INTEI.-KT  8141-028, 
April  1976. 


.^0 


INPUT  CARDS 


Card  Type 

Variable  Name  and  Dcscript ion 

lormat 

1 

Title  card  (colums  1-72) 

12A6 

HCALC(l)  - electron  emission  energies 
in  MeV.  Minimum  KCALC  is  0.001  MeV. 

Up  to  34  values.  Calculated  from  pho- 
ton energies  if  not  read  in. 

NAMIiLKST/SPbCIRA/ 

*> 

UVSPI'CT(l)  - energies  in  MeV  for  photon 
spectrum.  Up  to  50  values.  Omit  for 
blackbody  spectrum. 

NAMbLIST/SPLCTRA/ 

DNDfSP(I)  - relative  intensities 
(energy/energy  or  number/energy 
according  to  NSPPCT)  of  photon  spec- 
trum. Up  to  50  values.  Omit  for 
blackbody  spectrum. 

NAMbLlST/SPbCTRA/ 

KBT  = 0 (default)  if  an  arbitrary 
sjiectrum  is  input 
> 0 to  generate  an  incident 
blackbody  spectrum  with  a 
temperature  of  KBT  keV. 

N,\Mi;LIST/SPnCTRA/ 

9 

lOPT(l):  40  values 

I0PT(26)  = 1 to  punch  all  emission 
electron  energy  distri- 
butions on  cards  in  free- 
field  format  (default  = 0) 
I0PT(27)  = 1 to  put  photon  and  elec- 
tron spectra  on  file 

TAPK7  for  treatment  l)y 

PbOTAl.L  (default  = 1) 

NAME LI ST/ SPECTRA/ 

2 

NSPliCT  = 0 (default)  for  DNI)HSP(I) 
in  energy/energy 
= -1  for  I)NI)i:SP(l)  in 
number/energy 

NAMELIST/SPECTRA/ 

3 

INDHXX-QUlCKi;2  internal  table  number 
for  test  detector  element  for  dose 
calculation.  Must  specify  an  clement 
from  Table  4-1. 

WARNING:  DO  NOT  USE  ATOMIC  NUMBERS  HERE. 
SPECIFY  THE  APPROPRIATE  NUMBER  FROM  THE 
COLUMN  UBELED  "INTERNAL  NUMBER"  IN 

TABLE  4-1. 

15  (columns  1-5) 

^ NOTi: : 

( 

i 

Must  input  separate  NAME LI ST/ PLATE/  cards  for 
There  is  no  limit  on  the  number  of  i)lates. 

31 

each  plate. 

I 

I 


('arii  Type  Variable  Name  and  Descr i|)t  ion  l-ormat 

-1  NPCT  = 0 (default)  for  plate  elements  NAMF! LlST/l’LATi;/ 

in  fraction  by  number 
= 1 for  plate  elements  in  frac- 
tion by  weight 

4 II’RINT  = 0 (default)  no  cross-section  NAMld.IST/PIvMI./ 

|ir  i nt 

= 1 to  print  cross  sections 
from  input  tape 

4 KOPT  = 0 (default)  NAMlil-IST/PI^VlH/ 

= 1 to  print  detailed  information 
about  electron  production  by 
(ihotons 

4 NTYPF.  = 0 (default) 

= 1 to  calculate  bulk  current  at 
forward  edge  of  (ilate 

4 NKL(l):  atomic  numbers  of  elements  com-  N.AMliLIST/PLATF/ 

prising  plate.  See  Table  4-1.  Maximum 
number  = 8 per  plate. 

4 T ( 1 ) .'  respective  fractions  of  elements  NAMTl.lST/Pl.ATF./ 

(see  NPCT).  Default  = 1.  Maximum 
number  = 8 per  plate. 

4 lOPT(I):  40  values.  No  values  need  to  NAMFLlST/PL/VrF/ 

be  read  in  here  at  present. 

4 IFM  = 0 no  emission  from  plate  NAMFLIST/PLATF/ 

= 1 forward  emission  from  plate  (dcfaul t ) 

= 2 reverse  emission  from  plate 
= 3 forward  and  reverse  emission 
from  plate 

4 TAD:  thickness  of  plate  in  g/cm*"  NAMF.LlST/PlJkTF/ 

4 lIVOFtl:  angle  in  degrees  between  surface  NAMELlST/PlJVTF/ 

normal  of  plate  and  direction  of  photon 
beam.  Default  = 0. 

KARNlNd:  THIS  ANGLE  IS  NOT  INCLUDED  IN 
PHOTON  atti;nuation  calculation,  also, 

WHEN  NON-ZERO  VALUE  IS  USED,  RESULTING 
ELECTRON  YTEEDS  ARE  GIVEN  IN  UNITS  OF 
C/cm^  NORMAL  TO  MATERIAL  SURFACE;  PER 
cal/cm^  PERPENDICULAR  TO  PHOTON  BF.AM. 

4 QUIT  = 0 (default)  NAMELIST/PUTE/ 

= I signifies  last  plate 

DENSITY:  density  of  plate  (g/cm^) 


4 


N/XMELLST/PIJVTE/ 


i;i,i.Mt:NTS 

Table  4-1 

AVAlLABl.i;  IN  QU1CKT.2 

Atomic  Number  (Z) 

Internal  Number 

Name 

1 

1 

Hydrogen 

4 

Beryl  1 ium 

b 

3 

Carbon 

7 

4 

N i t rogcn 

8 

5 

0\ygcn 

9 

() 

I- 1 uor  i de 

12 

n 

Magnesium 

13 

8 

A 1 um i num 

14 

9 

Si  1 icon 

17 

10 

(dilorinc 

18 

11 

Argon 

20 

12 

Ca 1 c i um 

22 

13 

T i t an i um 

26 

1 1 

1 ron 

28 

15 

Nickel 

29 

16 

Copper 

32 

17 

Germanium 

41 

18 

N i ob i um 

47 

19 

S i Iver 

50 

20 

Tin 

. 60 

21 

Neodymium 

; 73 

22 

Tantalum 

79 

23 

Gold 

1 82 

24 

Lead 

92 

25 

Uranium 

S.  S,\MI'1J-;  PROBLKM  AM)  OU'll’DT  DliSCU  I I'T  I ON 


The  QUICKh2  out]Hit  is  best  illustrated  by  sample  jiroblems.  Samjile 
calculations  are  discussed  here,  i;iving  the  physical  and  cal  cu  1 at  i ona  1 
sij>n  i ficances  of  the  major  quantities. 

The  first  problem  involves  the  emission  j>enerated  when  a 15-keV 
hlackbody  jihoton  spectrum  is  incident  on  40  mils  of  M.  I'he  injiut  card 
images  for  this  jiroblem  arc  as  shown  in  Figure  5-1,  and  the  output  as  in 
Figure  5-2.  On  the  first  page  of  out[)ut  is  the  title,  followed  below  by 
a description  of  the  first  material,  or  "plate,"  since  the  problem  is  one- 
dimensional. This  description  gives  the  number  of  elements  in  the  mate- 
rial and  the  name,  fraction  by  number,  and  atomic  number  of  each  element. 
The  data  file  location  where  the  cross  sections  for  each  clement  are 
stored  is  given  (here,  "8"  indicates  that  the  infonnation  for  AH  is  stored 
on  file  8 of  the  data  tape),  and  the  thickness  of  the  material  is  jirinted. 

The  quantity  labeled  "input  flucnce"  is  simj)!)'  a convenience  output  of  the 

•) 

code.  It  gives  the  total  fluence  in  cal/cm"  of  the  photon  spectrum  as 

2 

input  to  the  code.  The  code  renormalizes  this  fluence  to  1 cal/cm*". 

An  important  point  regarding  QUlfKF2  output  should  be  understood  when 
interpreting  results  from  other  than  normally  incident  photon  beams.  Quan- 
tities given  in  terms  of  unit  area  differ  in  definition,  depending  on 
whether  they  pertain  to  photons  or  electrons.  Photon  densities  are  always 
per  unit  area  perpendicular  to  the  beam,  whereas  electron  densities  are 
always  per  unit  area  peri)endicular  to  the  surface.  Thus,  for  a non-normal 
incidence  problem,  the  correct  yield  is  obtained  by  simply  emitting  the 

number  of  electrons  per  unit  area  on  the  skewed  surface  obtained  directly 

2 

from  QUICKFZ  and,  of  course,  multiplying  by  the  flucnce  in  cal/cm  perpen- 
dicular to  the  photon  beam. 
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sample  problem  input  tor  15-kcV  blackbody  spectrum 
incident  on  a 40-mil  M plate  causing  both  forward  and 
reverse  electron  emission 
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Figure  5-2.  Sawple  problem  output  for  15-keV  blackbody  spectrum 
incident  on  40  mils  of  A£ 
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[When  reverse  emission  only  (emission  code  2)  or  forward 
and  reverse  emission  (emission  code  3)  is  asked  for  from 
a plate,  the  bulk  currents  just  inside  the  plate  face 
closest  to  the  source  are  printed  out.  Unfortunately, 
in  some  versions  of  the  code,  the  plate  index  number  for 
reverse  emission  is  low  by  I.  In  later  versions,  this 
error  has  been  corrected  and  the  legend  says  "Bulk  cur- 
rents in  layer  = at  entry  face."] 

MUL*  loXKlNtS  l^  s 1 {"ISSIU'.  coot  = 1 

t)UL«  COKXtN’S  »LLXC.  iMt 

M t bt'L"  LuxxlMS  PlxtCTIU'.t  N.^-uluOf-OB 

BtWfWSL  LuBxtM=  1 . 7SX|  ,;t  .OK 

HUL«  ruXXt\!s  MtXBk\(ilctL»t<  TCI  tHt  MJWMJL  PlxtrUtiS 
fOXoAXt/  COHXIM  I .',^,'(•1x1-0" 

BtViXSt  (UXX^^t  c 

[When  forward  emission  is  requested  from  a plate,  either 
by  itself  (emission  code  1)  or  with  reverse  current,  the 
bulk  currents  in  the  plate  near  the  face  farthest  from 
the  source  are  printed  out.  The  layer  index  is  correct. 
In  later  versions,  the  legend  says  "Bulk  currents  in 
layer  = at  exit  face."] 
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Figure  5-2  (cont . ) 
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In  tlio  next  section  of  out|)ut  , the  incident  photon  spectrum  normal- 

i;ed  to  a total  incident  fluence  of  1 cal/cnr  is  printed  out,  followed  by 

the  total  energy  content  of  the  spectrum,  the  average  |)hoton  energy,  and 

•> 

the  total  number  of  photons/cm“  corresjiond i ng  to  the  total  energy  con- 
tent. (Here,  the  total  numbei'  of  photons  is  1 calorie's  worth  of  jihotons 
with  the  spectral  energy  distribution  indicated.]  1 he  average  photon 
energy  is  taken  as  that  energv’,  hv,  which  gives  the  total  etierg.v  content 
when  multii'lied  by  the  total  number  of  incident  jihotons.  The  total  num- 
ber of  incident  jihotons  is  useful  in  that  results  from  other  codes 
exjiressed  in  e 1 ec  t rons/])hot  on  can  easily  be  converted  to  cou  1 ombs/ ca  1 or  i e 
by  multiplying  them  bs'  this  total  number  and  the  electronic  charge: 

Yield  in  C/cal  = yield  in  e/photon  x (photon/cal)  x (H/e)  . 

The  next  jiage  of  output  contains  elccti'on  emission  information.  The 
first  line  indicates  the  type  of  yield  there,  reverse),  the  material  num- 
ber from  which  the  electrons  are  being  emitted  (here,  1 indicates  the  first 
material),  and  the  yield  in  C/cal  incident  on  i)late  1.  This  yield  cor- 
responds to  the  "quantum  efficiency."  I'he  yields  in  C/cal  inciuent  on 
plate  .N  are  also  provided  in  case  the  user  is  interested  in  the  yield  pet- 
calorie  for  a filtered  spectrum. 

•Next,  the  breakdown  of  the  electron  emission  into  angular  and  energy 
distributions  is  given.  The  "emitted  electron  anguhir  distribution  inte- 
grtited  over  energy"  is  meant  to  provide  the  angular  distribution  only 
roughly.  The  "emitted  electron  energy  vli st r i but  ion  integrated  over  angle" 
should  be  a more  precise  function,  provided  the  energy  bins  are  close 
enough  together.  ibis  ])rintout  is  derived  from  the  "energy  distribution" 
described  in  Section  2.2.  The  column  KTIAI.  contains  the  number  of  elec- 
trons emitted  with  energies  greater  than  the  corresponding  entry  in  the 
column  h.NLRCY.  The  column  I’hK  BIN  is  the  total  number  of  electrons  in  a 
given  energy  bin  bounded  by  the  energies  in  the  column  BNliRGY  and  centered 
on  the  corresponding  energy  in  the  column  MIDI’Ol.NT.  NllMBlTt/MI'V  is  the 
number  of  elect rons/MeV  emitted  at  the  corresponding  energy  under  MIDPOINT. 
I’rinted  below  the  energy  distribution  is  the  average  electron  energy 
(calculated  in  the  same  way  as  the  average  photon  energy)  and  the  total 
emitted  charge.  The  latter  number  is  the  result  of  integrating  the  energy 
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1.1  i St  r ibut  ion  tlN/dli  over  energy.  It  should  agree  with  the  yield  at  the  top 
of  the  page  to  within  about  20  percent.  Agreement  is  generally  good 
except  for  back  emis.sion  produced  by  low-energy  photons  (about  5 to  10  keV'J. 
In  these  cases,  results  should  be  used  cautiously. 

^ fhe  next  two  pages  are  printer  plots  of  TOTAl.  versus  nNHRGY  and 
NtlMBliR/Ml.V  versus  MIOPOINT,  described  above.  These  jilots  are  exactly 
the  size  of  st.indard  S x 7i  log- log  iiajicr  and  may  be  xeroxed  or  used  with 
a light  table  to  save  needless  hours  of  hand-|)  1 ot t i ng . 

Ihe  cycle  repeats  again  with  a printout  of  the  photon  spectrum  emerg- 
ing from  the  first  material  (and  incident  on  the  second).  In  the  siunple 
problem,  the  low-energy  photons  have  been  absorbed  by  the  aluminum,  reduc- 
ing  the  total  energy  content  from  1.0  to  0.8J  cal/cm"',  raising  the  average 
photon  energy  from  40.5  to  50.8  keV,  and  reducing  the  total  number  of  pho- 
tons from  6.44  x 10^'  to  4.21  x 10^“  cm 

Dose  information  is  printed  at  this  point.  The  energy  depositions 
are  in  units  of  rads  of  the  particular  material  indicated  (1  rad  = 100 
‘•'I'R/r)-  Silicon  and  carbon  are  always  given,  and  a third  material  is  spe- 
cified in  the  injiuts  from  among  those  available  in  QUlCKi;2  (see  Table  4-1). 
Doses  at  the  fronts  and  backs  of  the  plates  are  given. 

The  forward  electron  emission  results  are  in  the  same  format  as  the 
reverse  ones,  but  the  spectrum  used  is  that  incident  on  the  second  mate- 
rial - i.e.,  the  spectrum  seen  by  all  material  within  a few  electron 
ranges  of  the  emitting  surface.  The  yield  from  plate  1 "per  cal  incident 
on  plate  2"  is  higher  than  the  yield  from  plate  1 "per  cal  incident  on 

plate  1"  since  the  latter  number  is  calculated  for  a total  spectral  energy 

2 

content  of  0.82  cal/cm  - i.e.,  the  energy  remaining  after  the  incident 
photons  have  been  filtered  through  40  mils  of  aluminum.  The  total  yields 
at  the  top  and  bottom  of  the  page  agree  to  4 percent,  so  in  this  sense, 
the  results  are  trustworthy. 

Bulk  current  information  is  also  given.  These  currents  are  in  units 
2 2 

of  electrons/cm  -sec,  photons/cm  -sec,  or  electrons/photon.  They  pertain 
to  regions  greater  than  the  largest  electron  range  away  from  any  interface. 
The  quantities  are  described  in  some  detail  in  Section  2. 


Finally,  for  extra  information,  there  are  two  |)rintinR  options.  If 
the  parameter  KOPT  is  set  to  1 (see  Section  J),  the  code  will  generate 
the  original  QI)1CKF.2  output.  If  IPRINT  is  set  to  1,  information  from  the 
data  tape  will  be  printed.  Information  on  these  printouts  is  available 
in  Reference  1. 
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(>.  Q1IICKF.2  COMI’UTFR  RF.QIJ  1 RFMli.NTS 


QU1CKF2  is  a FORTRAN  IV  computer  program  of  approximately  2000  cards 
in  length,  operational  on  a (’DC  7()00  computer.  No  machine  language  is 
employed.  The  code  requires  a tape  containing  the  tabulation  of  a con- 
siderable amount  of  data  for  material  properties  for  the  25  elements 
available.  This  file  is  about  .500,000  numbers  in  length. 

One  additional  file  is  also  required  for  plot  information.  Memory 
requirements  are  about  .54,000^^^  words  of  small  core  and  no  large  core. 

Typical  run  times  are  about  20  sec  CDC  7600  time  per  electron  emis- 
sion spectrum.  No  additional  computer  programs  are  required  by  QU1CKF2, 
but  the  PLOTALl.  code  (Ref.  11)  is  very  helpful  in  analyzing  results  of 
QIJ1CKF2  calculations. 

.A  present  idiosyncrasy  of  the  ()U1CKF2  version  described  here  is  that 
it  is  restricted  to  the  (!1)C  FTN  I'ORTRAN  extended  comjiiler,  version  1. 
Minor  programming  features  do  not  permit  its  use  on  the  more  up-to-date 
FTN4  compiler  available  currently. 
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